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Abstract 

The baryon self-energies are expressed in terms of the QCD condensates of the lowest dimension 

^— i in symmetric and asymmetric nuclear matter within the QCD sum-rule approach. The self-energies 

^^ are shown to satisfy the Gell-Mann-Okubo relations in the linear SU{?>) breaking approximation. 

bJQ The results are in qualitative agreement with those obtained by the standard nuclear physics 

■^ methods. 
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I. INTRODUCTION 

The study of the in-medium interactions of the octet of baryons is one of the hot topics 
in nuclear physics. While many years of theoretical and experimental investigation provided 
a very precise phenomenology of nucleon interactions, the hyperon properties in nuclear 
matter remain much less known. 

The experimental information on the hyperon in-medium interactions mainly comes 
from the hypernuclear physics. In the past years, an impressive experimental data on 
A-hypernuclei has been accumulated providing a potential depth Da ~ —30 MeV at the 
saturation density po = 0.17 fm~'^ [IJ. On the other hand, the unavailing search for S hy- 
pernuclear states [2j and the study of S^ atoms [3r5j show strong evidence for a repulsive 
nature of the S hyperon potential in nuclear matter. The S nuclear interactions seem to be 
attractive with the potential [/= ~ —18 MeV [6^. Finally, the hyperon-hyperon interac- 
tions were not really measured, there is just a handful of double A hypernuclear events [9]. 
Of course, the hypernuclear data is limited to the isospin-symmetric matter at the saturation 
density. 

However, the hyperon in-medium potentials are essential for the determination of the 
composition of neutron star matter p^fTT] . For example, in case of an attractive S hyperon 
potential, the S~ can appear even before the A hyperon in dense matter, but if S hyperon 
potential is repulsive then S hyperons are not populated at all. The hyperon population 
largely influences the mass-radius relation and maximum mass of neutron stars, the cooling 
of neutron stars, the stability with regard to the emission of gravitational waves and the 
possible early onset of the QCD phase transition in the neutron star cores [10]. Finally, the 
knowledge of hyperon in-medium properties is required to investigate an exciting possibility 
of strange hadronic matter [T2| stable against strong-interaction decays. 

There are several theoretical approaches to the problem of hyperon interactions in nuclear 
matter. The traditional method is based on the relativistic mean-fleld approximation (RMF) 
with effective meson-hyperon couplings flxed from the hypernuclear data and supplemented 
by the flavour SU{3) considerations P^ |T3]. A wide range of predictions was obtained in 
this approach, depending on the set of parameters chosen to describe the available data. 

Another commonly used method is based on the Brueckner-Hartree-Fock (BHF) approx- 
imation [T5HT7] with the soft-core hyperon-nucleon potentials [HI [19] extracted from the 



YN scattering experiments. In general, BHF approach successfully reproduces available 
hypernuclear data, however, the uncertainties in the scarce hyperon scattering data lead to 
large differences in the resulting YN in-medium potentials. 

One more method for the calculation of the A and S mass shifts is provided by the 
chiral perturbation theory (ChPT) [20]. However, this approach is limited to low densities 
up to 0.4po, and, even in this region, uncertainties are high since one has to determine 12 
low-energy constants from the fits to the scarce experimental data on the hyperon-nucleon 
scattering. 

Finally, there is a possibility to study the density dependence of the hyperon properties 
within the QCD sum rule approach which is based on the dispersion relations for the cor- 
relation functions of corresponding hadronic currents. Initially, QCD sum rule approach 
was developed to express vacuum characteristics of mesons in terms of expectation values of 
QCD operators known as "condensates" [21]. The sum rule approach was also applied to the 
calculation of nucleon properties in vacuum [22l [23] . Later, it was successfully extended to 
the studies of nucleon self-energies in nuclear matter [231 - I26] as well as to the calculation of 
the nucleon-nucleus scattering amplitude [27] and in-medium modifications of vector meson 
properties [2H]- An alternative version of the finite-density QCD sum rules, based on the 
Lehmann representation for the Green function, has been also developed [22] and applied 
to the calculation of A and E hyperon properties in symmetric nuclear matter [30ti33] . In 
contrast to other nuclear physics methods, the finite-density sum rule approach does not rely 
on phenomenological baryon parameters. The baryon in-medium properties are expressed 
in terms of the QCD condensates, which can be either calculated in a model-independent 
way or related to observables. 

In this paper, the finite-density QCD sum rule approach [26] is extended to the calculation 
of the baryon octet properties in nuclear matter. The baryon self-energies are expressed in 
terms of the lowest dimension quark and gluon condensates taken in the gas approximation. 
It is known that the SU{'i) symmetry breaking in the baryon octet in-medium properties is 
caused both by the nonvanishing strange quark mass iris and the SU{'i) asymmetric quark 
composition of the medium itself. We will show in this paper that the baryon self-energies 
satisfy the Gell-Mann-Okubo relations in the linear SU{?)) breaking approximation. The 
numerical results are obtained for the symmetric and asymmetric nuclear matter and studied 
as functions of the scalar quark condensate. Compared to papers [301433] . we concentrate on 



the effects of the broken SU{3) symmetry in the baryon octet and extend the hyperon sum 
rules to the case of asymmetric nuclear matter providing a framework for the calculation of 
the neutron star equation of state within the QCD sum rules approach. The properties of 
the S hyperon in nuclear matter are considered for the first time. 

The structure of the paper is as follows. In Sec. [TT| the QCD sum rules for the baryon octet 



in vacuum are reviewed. In Sec. |III| the finite-density QCD sum rules are generalized to the 
case of the matter consisting of an arbitrary mixture of baryons, the approximate solutions 
are obtained and the numerical results in the symmetric and asymmetric nuclear matter 



are discussed. In Sec. IV, the results are compared to experimental data and predictions of 



other approaches. A summary is provided in Sec. IVl 



II. BARYON OCTET IN VACUUM 



The QCD sum rule approach in vacuum is based on the dispersion relation for the cor- 
relation function 

T^Boiq') = ^Jd^xe'^-'{0\rjBix)MO) |0), (1) 

where T denotes the time-ordered product and Jb is the local three-quark current with 
quantum numbers of the baryon in interest. The usual choice for the proton current is [22l 



jp = eabciu^C'yf,Ub)'y5'l''dc , (2) 

where a, b and c are the color indices, T denotes the transpose and C is the the charge 
conjugation matrix. The currents for other members of the baryon octet may be obtained 
by the flavour SU{3) transformation of the proton current [35] . 

The general idea of the QCD sum rules is to approach the bound state problem in QCD 
from the asymptotic freedom side. At large negative g^, the correlation function IIbo{q^) 
is approximated by the power series in g^^ known as the Operator Product Expansion 
(OPE) [36]. On the other hand, the imaginary part of IlBo{q^) at g^ > can be described 
in terms of the observable hadrons. Two momentum regions are connected in the dispersion 
relation for the function IlBo{q^): 

^Boici) = -J^^^—^dk. (3) 



At large negtive g^, the left-hand side of this equation is approximated by several lowest 
terms of OPE n^ol?^) ~ n^Q^(g^) with the coefficients containing the expectation values 
of the local quark and gluon field operators. 

The phenomenological right-hand side of equation (|3| is usually considered in the "pole 
-|- continuum" model where the lowest lying pole in iTaUsoik^), corresponding to the baryon 
in interest, is separated from the higher k"^ singularities approximated by the continuum [21] : 

ImUsoie) = Xl6{e - ml) + ^O^e - W^|)Ang™(A;2) (4) 

where A^ and tub are the baryon residue and mass while Wb represents the effective con- 
tinuum threshold value. An§Q^(A;^) denotes a discontinuity caused by logarithmic terms 
of the perturbative expansion of the correlation function. Thus the dispersion relation ([3]) 
takes the form: 

tOPE/„2n _ ^B , ^ f^ ^n^Q [k ) 2 



rriB — q 2m jw% k^ — q^ 

The perturbative expansion at the left-hand side becomes increasingly valid for large 
|g^| = Q"^, while the "pole + continuum" assumption becomes more accurate when |g^| 
decreases. Usually, a certain intermediate region of g^ values is considered where approxi- 
mations on both sides of equation ^ are believed to be valid. To improve the overlap of 
two approximations, the Borel transform is usually applied to both sides of equation ^ 
converting a function of Q^ into the function of the Borel mass M^ |21j . 

The correlation function IlBo{q^) can be decomposed into two structures: 

UBoiq') = qU%{q') + mUq') , (6) 

where / represents the identity matrix and q = q^'j'^. The Borel transformed dispersion 
relations for the structures n^Q(g^) and n^Q(g^) are known as QCD sum rules in vacuum 

AUUk'' 



KoiM') - 1"^ ""'^ ^ e-''/^''de = \\e-^y^' (7) 

nio(M^) - r ^Hk(^e-^'VM^^^2 ^ ^^A|e-B/M^ (8) 

Jw% 2-111 

It is convenient to write these equations in a compact form: 

C'b,{M\ Wl) = nUM', TUB, A|), l = q,I (9) 



where C''^q{M'^, W^) represent the left-hand sides of the sum rule equations ([TJ-JS]) multiplied 
by the factor 327r^ while the right-hand sides are expressed as 






(10) 
(11) 



with A^ 



327r^A|. 



The sum rules for baryons are usually considered in a certain interval of the Borel mass 
M^ where the contribution of the continuum and higher OPE corrections are found to be 
relatively small [37]: 

0.8 GeV^ < M^ < 1.4 GeV^ (12) 



The unknown values rriB, Xb and W^ are obtained by minimization of the function 



x\mB,Xl,Wl) 



J l=q,I 



'CUM,)-KoiM, 



(13) 



^U(M,) ; ' 

which insures the most accurate approximation of C^Boi^"^^ ^b) by T^mi^'^y ""^-B, A^) with 



Mj being a set of points evenly spaced within the fiducial interval (12). There is also an 



alternative approach based on the logarithmic measure [371 EH] • 

Initially, the QCD sum rule approach was applied to the calculation of the baryon masses 
in [221 1231 Ell EH]- in the nucleon case, the following expressions were obtained for the 
left-hand sides: 



C%,{M',W'^) 



cUm',w^ 



B'ia + B-jah + Bo,c? , 



(14) 
(15) 



where /Xq = 0.8 GeV^, while for quark and gluon condensates the traditional notations are 
used: 



a = -(27r)2(0|gg|0), 6 = (27r)Vo 



ftc 



G' 



TX 







(16) 



The numerical value for the quark condensate (0|gg|0) = —(0.24 GeV)'^ was obtained from 
the well-known Gell-Mann-Oakes-Renner relation [40j while the value of the gluon con- 
densate (0 ^G"^ 0\ = (0.33 GeV)"^ was extracted from the analysis of leptonic decays of 
p and mesons and supported by the QCD sum rule analysis of the charmonium spec- 
trum [H]. Note that isotopic invariance is assumed for the light quark condensates with q 



denoting the light quark field, while higher-dimensional quark condensates are considered 

in the factorization approximation: {0\qqqq\0) = ((0|gg|0))^ pi] . 

The coefficients A„ and -B„ are functions of M^ and W^ with the subscript denoting the 

dimension of the corresponding condensate: 

M^E2 _ M'^Eo , 4.,/n . 1 



An 



A, 



Bs = 2M^Ei, 



Br = — 



1 



Bo 



3 

272a, 



A, 



37^^-2^2/27' 



9' "" 8l7rM2LV9' 
These expressions depend on the continuum threshold via the functions En defined as: 



E„ 



1-e-^E 



k=0 



X 



X 



M2 



(17) 



(18) 



Finally, the function L = L[M'^) accounts for the leading logarithmic corrections |21j : 

Here A = Aqcd = 0.15 GeV, while v = 0.5 GeV is the OPE normalization point. 

The expression for the nucleon mass as function of Borel mass M and continuum threshold 

Wn directly follows from ^, (14) and (15): 

//i^2 T;r/2N B^a + Bjab -\- Bga 

mN{M^,Wtr) = — — 5"^ (19) 

^ ^^ Ao + A^b + Aea^ + Asfx^a^ ^ ' 

Note, that this expression was obtained in the chiral SU{2) limit (?7i„ = rrid = and 
(OImmIO) = {0\dd\0)). In the case of hyperons, the strange quark mass rris and the difference 
in values of strange an light quark condensates become important. The value of the strange 
quark mass rus is about 150 MeV with the uncertainty of 20% [12]. The deviation of 
the strange quark condensate from the light quark condensate is usually described by the 
parameter 7: 



7 



{0\ss\0) 



-1. 



(20) 



{0\uu\0) 
The value 7 = —0.2 is usually accepted [3S1 EHl EH]- Following general QCD sum rule 

technics, one may obtain the mass formulas for hyperons with the condensates accounted 

up to dimension 9 [39lH3lHi]: 



rriA 



ms 



rri'z 



(53 + Bjb) (1 



fia 



59(1 + 7)«' 



Sn - S.b - 4Sa(l 



IW 



Ao + A,b + {A, + Asfil) (1 + f ) a2 + [1^3(1 - 87) - ^5(1 - l)f^l 

{B3 + Brb){l + 7)a + ^9(1 + 7)a' + ('^0 - '^4^ + ^Sea^) m, 
Ao + Aib + (Ag + A^^il)a^ - (^3 + 55/ig)(l + -f)ams 
{Bs + Bjb)a + ^9(1 + 7)'a3 + 3^6(1 + 7)q'^. 
Ao + A^b + (Ag + A8/ig)(l + -ffa? - 2^5(1 + 7)/igam, 



m.o 



am. 



(21) 

(22) 
(23) 



where Si are coefficients in the additional terms hnear in m^: 

2M^E',. „ M^Eo ^4 ^ 2M'^Eo 



So 



Sa 



S, 



S. 



S, 



1 



(24) 



^8/9 ' ^^ 4L8/9 ' "^ 3' '"■' L'^/^ ' """^ 3L26/27- 
Considering the chiral SU{3) hmit {'~f,ms — ?■ 0), one can check that the expressions for 



hyperon masses (21-23) reduce to the the nucleon mass formula (19). Moreover, GelL 



Mann-Okubo mass relation 



2(m7v + rn^) = m^ + 3mA, 



(25) 



which is a direct consequence of the symmetry breaking in the (3, 3) term of the SU{3) 
Hamiltonian, is automatically satisfied up to the terms linear in m^ and 7 |39j. However, 
one has to keep in mind that this property is valid only in the limit of equal continuum 
threshold values W^ = Wl = Wl = Wl 

The parameters m^, A^ and W^, which minimize the function (13), are shown in Table ^ 
The obtained baryon masses agree with the experimental values within 10% accuracy. The 
systematic underestimates can be attributed to the fact that radiative corrections |13] or 
non-perturbative effects due to instantons [l6l HT] may be important. 



B 


m^^P, GeV 


TUB, GeV 


A|, GeV^ 


W^, GeV^ 


N 


0.940 


0.934 


1.897 


2.119 


A 


1.116 


1.103 


3.189 


3.069 


S 


1.193 


1.104 


3.066 


3.157 


S 


1.314 


1.207 


4.069 


3.729 



TABLE I. Values of m^, A^ and W'^ from the mininiization procedure. 

Variation of the strange quark mass and the parameter 7 results in significant changes 
of the hyperon masses and their relative order. Therefore, it is instructive to consider 
approximate expressions for the hyperon masses as functions of the involved parameters. 
Multiplying both sides of the QCD sum rules by exp(m^/M^) one may notice that the 
functions C''qq{M'^, l^^)e™sM^ should not depend on the Borel mass. A simple check ensures 
that the leading OPE terms, multiplied by the factor exp(m^/M^), can also be approximated 



by the constants in the range (12) within 10% accuracy. Referring to this observation, let 

us introduce the notation: 

X„(M2,iy|)exp(m|/M2) 



Xrr{mB,W^,X%) 



\2 



(26) 



where X„ stands for the functions An, Bn or S'„, while the overhne denotes averaging over 



the Borel mass range (12). The numerical values for A^, Bn and Sn were calculated at values 



of niBi W^ and A^ from the minimization procedure and provided in Table M Then the 



mass formulas (19), (21), (22), (23) can be expressed in terms of the averaged values An, 
Bn and Sn instead of M^-dependent functions. Such expressions reproduce the values of 



the baryon masses, fitted via x^ (13), with an accuracy of the order of 0.5%. For example. 



in the case of nucleon and H hyperon, one can write: 

1.61a -0.10a6 + 0.43a3 



ruN 



m'= 



0.38 + 0.226 + 1.97a2 - O.ST/iga^ ' 

1.91a - 0.08a6 + 0.36(1 + 7) V + 3.93(1 + -i)a^ms 
0.65 + 0.206 + (1.61 - 0.31/ig)(l + 7)^ - 0.42(1 + 7)/igam, 



(27) 
(28) 



where all values are expressed in powers of GeV. To our knowledge, such expressions for 
the baryon masses are considered for the first time. They reveal relative contributions of 
different OPE terms providing a convenient way to study the dependence on the condensate 
values. Similarly, using the expressions for the hyperon masses, it is easy to show that 
reasonable variations of irts and 7 values do not allow tuning all hyperon masses to their 
experimental values. Therefore, we will apply conventional values of irts = 150 MeV and 
7 = —0.2 in the finite-density QCD sum rule analysis. 



B 


^0 


A, 


Aq 


^8 


^3 


B7 


B9 


So 


^3 


S, 


S5 


^6 


N 


0.38 


0.22 


1.97 


-0.37 


1.61 


-0.10 


0.43 


- 


- 


- 


- 


- 


A 


0.54 


0.20 


1.63 


-0.31 


1.76 


-0.08 


0.36 


0.87 


1.61 


0.16 


0.21 


1.32 


S 


0.58 


0.21 


1.70 


-0.32 


1.87 


-0.09 


0.37 


0.94 


1.69 


0.17 


0.22 


1.38 


^ 


0.65 


0.20 


1.61 


-0.31 


1.91 


-0.08 


0.36 


- 


- 


- 


0.21 


1.31 



TABLE II. Values A„, Bn and Sn averaged over the Borel mass range according to equation (26) 
in powers of GeV. 



III. BARYON OCTET IN NUCLEAR MATTER 



In this section, we will develop the framework for the calculation of the baryon octet 
parameters in nuclear matter following the finite-density sum rule approach reviewed in |26] . 



A. QCD sum rules in nuclear matter 

The propagation of a system with four-momentum q in nuclear matter is described by 
the correlation function: 

I^Bmiq^) = 1 1 d^xe''^-^M\ TjBixyjBiO) |M), (29) 

where \M) is the ground state of nuclear matter. Considering the nuclear matter as a system 
of A nucleons with momenta Pi, one may introduce the vector 

P = ^ (30) 

which turns to p ~ {m, 0) in the rest frame of the matter, where m is the nucleon mass. 
The spectrum of the function YlBmiq'^) appears to be much more complicated than that of 
11^0(0'^)) however by fixing the value oi s = {p + g)^, one may separate the singularities 
connected with the matter itself from those connected with the baryon in the matter [25]- 
[27] . In this paper, the effects of the nucleon Fermi motion on the baryon properties will be 
neglected, therefore the threshold value oi s = {m + nisY will be used in the calculations. 
The general form of the polarization function in the nuclear matter can be presented as: 



IiBm{q) = qTi%m{q\ s) + m's„M\ s) + ^^^^(g^ s). (31) 



The in-medium QCD sum rules are then derived as the Borel-transformed dispersion rela- 
tions for the components n^„(g^, s)\ 

T^'Bm{<i\s) = -l^^^§^^^^dk\ ^ = q,I,p. (32) 

It was shown that the spectrum of the function nBm('?^,s) can be described by the 
"pole-|-continuum" model similar to the vacuum case at least until the terms of the order 
p^ are included in the OPE [26]. One may consider a general expression for the propagator 
of the baryon B in nuclear matter: 

Gb' = {G'b)-' - Sb, (33) 

where G% = {q — rnB)~^ is the free baryon propagator and S^ is a general expression for 
the baryon self-energy in nuclear matter: 

SB = gS^ + ^S^ + /Si. (34) 

m 

10 



Inverting G^ ^ we find for tlie in-medium baryon propagator: 









wliere Zb = [(1 — Sg)(l + S^/m)] \ wliile S^ and m^ correspond to tlie vector self-energy 
and the effective mass in nuclear physics: 

S^ = ^, -^ = ^^^^. (36) 

For the new position of the baryon pole niBm we find: 

^1- = ^ f, ^y/L ^ ■ (37) 



Following definitions, accepted in nuclear physics, it is also convenient to introduce the scalar 
self-energy E;| = m*^ — itlb and the non-relativistic baryon potential Ub = S)^ + S;|. 
The Borel transformed sum rule equations take the form (/ = g, /,p): 

C'Bm{M\WlJ=n'B^{M') (38) 

with the phenomenological right-hand side: 

^U(M^) = e'A|^e--/*^^ (39) 

where A^^ = 327r''ZBA^^ is the effective value of the residue for the baryon B in nuclear 



matter. The values of ^ are determined from equation (35): 



e = l, e = -S^, e = m*B. (40) 



The left-hand sides of the sum rule equations (38) are calculated in the OPE approach: 



^'Bn.{M\ WlJ = 32vr^ {n^Bm{M\ s) - Q ^^!n|^^e-^V^^^rfp) . (41) 

Similar to the vacuum case, we can express baryon effective masses and vector self-energies 
via the left-hand sides of the sum rule equations: 

r-MMM^lJ = #4^^^^, (42) 



C%^{M\W^bJ 



^Vf%f2 TTr2 \ ^Bm\^ i^Bm) 



^UM\ WU = - rr)^^JT\ - (43) 

'-BmV^^ ■> ^^ Br. 



imJ 



For the calculation of baryon in-medium properties, we will consider only leading OPE 
terms which density dependence is briefly reviewed in the next subsection. 

11 



B. Condensates in nuclear matter 

In this subsection, we consider the condensates in the nuclear matter of density p con- 
sisting of an arbitrary mixture of the baryon octet members B = p,n, A, S'^, S°, S~, S°, S~ 
with concentrations Cb- 

The lowest order of OPE in medium can be presented in terms of the vector and scalar 
quark condensates. The vector quark condensate is defined as: 

vUp) ^ {M\qa,q^\M), (44) 

where qt stands for u, d or s quarks. In the rest frame of nuclear matter, the vector con- 
densates take the form f^(p) = Vq{p)6^o, where the functions Vq{p) are linear in the nuclear 
matter density p: 

^o(p) = ^iP' ^i = J2^bCb (45) 

B 

with n^ = {B\qi'jQqi\B) denoting the number of valence quarks of flavour i in baryon B. 
For the ordinary nuclear matter consisting of protons and neutrons only, it is convenient to 
define isospin symmetric and asymmetric combinations 



V =Vu + Vd= {p\u-fou + d-fod\p) = 3, (46) 

v~ =Vu-Vd = {p\u'joU - d'jod\p) = 1. (47) 

While due to the vector current conservation the vector condensates are exactly linear in p, 
the scalar quark condensates K^(p) = {M\qiqi\M) are more complicated functions of density. 
However, in the gas approximation, they can also be expressed by the linear functions of p: 

f^lnip) ~ /^O + f^iP^ f^i = Yl f^B^B, (48) 

B 

where Kq = {0\qiqi\0) while k^ denote baryon matrix elements: 

4 = {B\Mi\B). (49) 



Following [49], we introduce isospin symmetric and asymmetric combinations of the light 
quark expectation values: 

K = k'^ + k'1,= {p\uu + dd\p), (50) 

C = Kp - Kp = {p\uu - dd\p) . (51) 

12 



The expectation value k is directly related to the pion-nucleon sigma term ct^at [50] : 

K = ^^^^ (52) 

with rriu ~ 4 MeV and rrid ~ 7 MeV denoting the current masses of the light quarks. The a 
term can be extracted in several ways, i. e. from the subthreshold extrapolation of the vrA^ 
scattering amplitude, however there is a large discrepancy between the results (see [25] for 
references). Assuming the conventional value of cr ~ 45 MeV [SI], one would obtain k ^ 8, 
while with the latest results of cr ~ 60 MeV, the value of k ?^ 11 is preferred [26] . 

In contrast, the expectation value of ( is not restricted by experimental data, therefore 
some model assumptions on the quark structure of the nucleon are required. If the nucleon 
is treated as a system of valence quarks and isospin-symmetric sea of quark-antiquark pairs, 
the expectation value ( is determined by the contribution of the valence quarks. We will use 
the value oi ( = 0.54 obtained in the perturbative chiral quark model (PCQM) [52] which 
was used in |48j for the calculation of in-medium four-quark condensates. 

The strange quark expectation value Kp = {p\ss\p) is usually parameterized in terms of 
the strange quark content y: 

{p\uu + dd\p) K 

The parameter y is strongly correlated with the value of the a term. The value cr ~ 60 MeV 
corresponds to the large y ~ 0.35, while the conventional value of a ~ 45 MeV is consistent 
with the smaller strange quark content y ~ 0.2. In PCQM approach [53j, one gets y = 0.08 
in support of the smaller strange quark content. We will use y = 0.08 as a default value in 
this paper since the values y = 0.2 or 0.35 would correspond to too many strange quarks 
in the nucleon in contradiction with the naive non-relativistic quark model. We will also 
study the sensitivity of the baryon self-energies with respect to the strange quark content 
parameter. 

The scalar quark condensate was also considered beyond the gas approximation in the 
framework of the meson-exchange model of nucleon-nucleon interactions [211 125] • It was 
shown that the nonlinear contribution to the scalar condensate may be responsible for the 
saturation mechanism. However, the nonlinear terms appear to be small compared to the 
linear term up to the saturation density, therefore they will be ignored in the framework of 
this paper. 

13 



The gluon condensate in nuclear matter can be also considered in the gas approximation: 

a.. 



9r. 



M 



'-a- 



vr 



M 



90 + 9P 



(54) 



where 5^0 
element: 



a^Q-. 



0) is the vacuum expectation value while g is the nucleon matrix 



9 



N 



"s^2 



'G' 



n 



N 



(55) 



The value of g was calculated in [M] by averaging the trace of the QCD energy-momentum 
tensor. In the chiral SU{3) limit, one gets g = ~^rn which is sufficient for our analysis 
referring to the small contribution of the gluon condensate [26] . 

As for the scalar quark expectation values {H\qiqi\H) in hyperons H = A,S, S, they 
cannot be directly related to observables and some model assumptions are necessary. One 
option is to apply Hellmann-Feynman theorem to the QCD Hamiltonian density and relate 
the scalar quark expectation value to the derivative dmH/dniq^ 

dniH 



mq,{H\qq\H) 



m. 



' dm. 



(56) 



li 



The functions niHifng.) and corresponding derivatives can be deduced from ChPT [32] or 
vacuum QCD sum rules. Note, however, that numerical results of this paper will be limited 
to the case of the nonstrange nuclear matter consisting of protons and neutrons only, thus 
we will not need the values {H\qq\H) in the present calculations. Nevertheless, our approach 
can be easily extended to the case of an arbitrary mixture of baryons. 

C. Sum rules in the gas approximation 



Following [26], we can express the left-hand sides (41) as a sum of vacuum expressions 
and terms linear in density p: 



^BmiM', WlJ = CUM', WlJ + X'j,{M\ WlJp, 



(57) 



where / denotes the structures q, I and p. Note, that £^q(M^,PV^^) = while the vac- 
uum expressions £^g are calculated at density-dependent continuum thresholds WBm- The 
functions X^ can be expressed in terms of quark and gluon expectation values Vi, Ki and g, 
considered in the previous subsection: 



Xi 



+ ^B Z^ 0,vB''^i + ''TT's^B Z^ O.KB'^i 



(58) 
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Xi = B^^Yl ^B^i + ^sBl E K^Vi. (59) 

i i 

X'b = PBJ:PlBV^. (60) 

i 

where the Borel transformed OPE coefficients A^, A^, A'"^, B'^, B"^ and Pg read: 

TiHl^ S7r ^s-m^)M'Eo-M^E, 

A^ = 47r2M2Eo, i?l^ = 47r2^ ^—^^ 1, (62) 

5^ = -4.VE„ P]^ = -^^, (63) 

These expressions depend on the Borel mass M^ and on W^^ via the functions Eo{W^^/M'^) 
and £'i(iy^^/M^) which are taken at density-dependent continuum threshold values. Note 



that the coefficients (61-63) also depend on s which was fixed in the dispersion relation 
for the correlation function. Analogous expressions were obtained in the finite-density sum 
rule approach based on the Lehmann representation for the Green function [301 EI] where 
the dispersion relation was written in the go complex plane at fixed three-momentum q. In 



the latter case, the Borel transformed OPE coefficients (61-63) appear to depend on the 
momentum q instead of the s invariant. 



The coefficients a^^, a^^, 6^^, 6^^ and p^^ are shown in Table III They depend on the 
baryon isospin projection I^b responsible for the splitting of baryon self-energies in the {n,p) 
and (S~, S°) isospin doublets and the (S^, S'^, S+) triplet. However, in the isospin symmetric 
nuclear matter, the terms, proportional to I^b, cancel out, and the isospin symmetry for the 
baryon self-energies is restored. 

Let us study the symmetry properties of the obtained expressions. It is easy to check that 
in the chiral SU{3) limit {rus — > 0), the SU{3) symmetry remains broken due to different 
coefficients aj,^, 6^.^ and pl^ accompanying light and strange quark condensates for different 
baryon species. However, the SU{3) symmetry is restored in the SU{3) symmetric matter 
with equal scalar and vector quark compositions {vu = Vd = Vs and Hu = i^d = i^s), providing 
degenerate functions X^ and hence equal effective masses (42) and vector self-energies (43) 



for the baryon octet. Besides, in the isospin symmetric matter, the functions X^ satisfy 
relations similar to the Gell-Mann-Okubo mass formula: 

2(X^ + XL) = X^ + 3X1 (64) 
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s 


^ 


<B 


1 


5 
6 


5(1 + ^3e) 


^(1 + 2/3h) 


<B 


1 
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6 


2(1 -^3S) 


2(1-2/3=) 


<B 





2 
6 


1 


1 


<B 





4 
3 








O'kB 





4 
3 








<B 
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1 





Kb 





1 
3 


-1-/3S 





Kb 





1 
3 


-1 + /3s 





Kb 





2 
3 


2 





Kb 


1 - 2/37V 


4 
3 





1 + 2/3= 


Kb 


1 + 2/37V 


4 
3 





1 - 2/3= 


Kb 





2 
~3 


2 





PvB 


1 + IhN 


11 
24 


i(l + /3E) 


|(1 + 2/3h) 


PvB 


1 - IhN 
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24 


1(1 -^3S) 


1(1-2/3=) 


pIb 





26 
24 


2 

8 


14 

8 



TABLE III. The coefficients a^^, a^^, 6J,^, h^^ and pj,^ in the expressions (58-60) 



which are vahd in the hmit of equal effective continuum thresholds for the terms, proportional 
to the gluon and light quark condensates, and for the terms linear in tti^, k^ or Vs- Under 
this approximation, the same relations also hold for the effective baryon masses and vector 
self-energies: 



2(m^ + mV) = m^ + 3m^ 



A' 



2(s](; + s^) = s^+3sr. 



(65) 
(66) 



Of course, these relations should be valid for any model based on SU{3) symmetry breaking 
hypothesis. 

Exact solutions for the baryon effective masses and vector self-energies, as well as in- 
medium effective thresholds W^m ^^^1 residues A^^, could be found by minimization of 



the function similar to (13). Note, however, that exact solutions of the nucleon sum rules, 
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accounting for four-quark condensates and nonlocalities of the lowest dimension conden- 
sates, resulted only in slight changes of the continuum threshold W^^ up to the saturation 
density [IB]. Therefore we can safely consider approximate solutions described in the next 
subsection. 



D. Approximate solution 



In this subsection, we consider an approximate solution for the sum rule equations (42) 
and (43) by replacing the in-medium continuum threshold W^^ by its vacuum value W^ [26] . 



r2 ,,.2^ /:iJo(M^ly|) + x^(MMy|)p 



We can divide both numerators and denominators by £^Q(M^,iy^) and express m*^ and 
S^ in the form which is generally accepted in the nuclear physics: 






s^(M^ wi) = - /^^,,,;;g;,, , (70) 



where 

Applying the vacuum sum rule, one gets: 



2^ 



^B 

Similar to the vacuum case, one can average the functions J-"^ over the Borel mass range at 
the vacuum continuum threshold and express them in terms of the average values for the 



functions Ag, A^, A^, B^, B^, P^: 



Tl = lA%g + AIY1 <BV^ + m,Al ^ a\j,KA p, (73) 

^i = (bb E KBf^^ + ^sBl E KbV^ P, (74) 

n=(nT.I^.BV^P. (75) 



I 
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where the meaning of the overhne was defined in (26). The values Ag, Ay, A^, By, B,^, Py, 



averaged at vacuum continuum thresholds for A^, A, S and S baryons, are shown in Table IV 



Then the effective baryon masses and the vector self-energies (69-70) can be expressed in 
terms of average values „F^ and m^: 



— * 
m 



7l 
B 



B 



i+n ' 



'V 



i+n' 



(76) 



B ^ ^ -> B 

In the next subsection, we will study these approximate solutions in the case of nonstrange 
nuclear matter. 
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A% 
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Al 
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-65 


-74 


-A% 
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39 


41 


40 


B% 


-32 


-35 


-37 


-38 


^B 


58 


68 


80 


90 


^B 


-69 


-75 


-80 


-82 



TABLE IV. The values Ag, Ay, A,^, By, B,^, Py, averaged at vacuum continuum thresholds (in 
powers of GeV). The numbers are round off to integer values. 



E. Symmetric nuclear matter 

Let us consider the nuclear matter, composed of protons and neutrons with baryonic 
concentrations Cp and c„ respectively. We can also define the isospin asymmetry parameter 
P = Cn — Cp, which is equal to 1 in the pure neutron matter. Then the functions X^ can 
be expressed in terms of the nucleon matrix elements v, v~ , k, (, and the strange quark 
content y: 



Xb = ^b9 + ^ij(a;!B^ + «,;B^ ^3b/3) + msABia^Bf^ + (^IbI^v), 
Xi = B%{htj,K + h-BChB^ + hlB^^y) + m,Bl{htj,v + Kj^v-hB^), 

X^B=PBiptBV+P^BV-hBP) 

h^Bi (^^Bi K.B ^^d PyB summarized in Table [v} 
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(77) 
(78) 
(79) 



with coefficients a^^, 6^^, a^^. 
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1 

8 


PvB 


3 

2 





7 
8 


1 
4 



TABLE V. Coefficients a^^, b^^, a+^, 6^^, 



^Ib^Kb andp^^. 



In the symmetric nuclear matter, all /3-dependent terms vanish, and the approximate 



solutions ( 76 ) can be written in a simple form: 



^* mg + [B^bKKb + KbV) + msBlhtBv]p 

^ 1 + {A's9 + ^Bfl+s + msA%K{atB + aLl/)]p' 
nv -PbP^b^P 



^B 



(80) 
(81) 



1 + [A%g + A^a^s + m,A%K{atB + alBy)]p' 
The obtained density dependence for the baryon effective masses and the vector self- 
energies is shown in Fig. [T]for the default values k = 8 and y = 0.08. One can observe that 
the ratios m*^/mB for A^, A and S follow almost identical trend and are around 0.8 at the 
saturation density. This coincidence is related to the fact that the main contribution to the 
effective masses comes from the term proportional to the coefficient h^^ which equals to 1 
for N and H and | for A (see Table M). In case of A, there are two additional terms in the 



numerator of equation (80), proportional to 6^^ and 6^^, which compensate the difference 



of K^K = 3 from 6^jv = K-e = 1- O^ ^^e other hand, the coefficient K^b ^^ ^ is equal to 
and the main contribution to the density dependence of the S effective mass comes from 
the term proportional to the strange quark content y and the strange quark mass mg. The 
interplay of these terms results in a positive slope of the S effective mass. 
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The slope of the vector self-energy E^, shown in Fig. ol is basically determined by the 
coefficient p^^ in the numerator of equation (81). According to Table M, the vector self- 
energies for A^, A, E and H should approximately scale as 1 : || : | : | in contract to the 
prediction 1 : | : | : | of the naive quark model. Different continuum thresholds result in 
different average values for the functions P^ which effects the deviation from the scaling 
: ^ : ^. The ratio of the nucleon vector self-energy to the vacuum nucleon mass is 



1 : 



11 

24 



about 36% at saturation density which is in agreement with results obtained in [39] for the 
lowest-order condensates. 
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FIG. 1. The density dependence of the baryon effective mass m*^ (a) and the vector seff-energy 
S}j (b) in the symmetric nuclear matter. 

It is also instructive to study the dependence of the baryon effective masses and vector 
self-energies on the expectation value n and the strange quark content y which appear 



to be somewhat ambiguous according to the discussion of section IIIB As for the vector 



self-energy, its dependence on the scalar quark expectation values is marginal according to 



equation (81), since it comes only from the term proportional to nis in the denominator. 
Moreover, in case of A^ and S baryons, this dependence completely vanishes since a^^ = 



''kN 



a" = 



0. Numerical analysis shows that any reasonable variations of k, and y 



do not change the values of the A and S vector self-energies within the accuracy of 0.5%. 

On the other hand, the parameters k, and y play important roles in the calculation 
of the baryon effective mass since they appear in the leading terms in the numerator of 



equation (80). Variation of the baryon effective masses versus k at the saturation density is 



illustrated in Fig. |2} One can observe that the effective masses for A^, A and S baryons drop 
down from 0.8mB to ~ O.GSm^ when k is varied from the conventional value k = 8 to the 
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value K = 11, which is favoured by the recent results on the nucleon a term. On the other 
hand, S effective mass has only marginal dependence on k since b'^j^ = 0. 
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FIG. 2. Baryon effective masses as functions of k at the saturation density. 

Finally, the variation of the baryon effective masses versus the strange quark content y 
is illustrated in Fig. [3j In case of A^ and S baryons, there is no dependence on y since 



VB 



bl 



kB 



for 5 = A^, H. The effective mass for the A hyperon changes by about 5% 



when the strange quark content y varies from to the somewhat extreme value of 0.35. 
On the other hand, the effective mass for the S hyperon dramatically depends on y since 
the strange scalar quark expectation value appears in the leading term due to vanishing 
contributions of the light scalar quark expectation values. Numerically, when y increases 
from to 0.35, the ratio m^jvaY. drops down from 1.1 to 0.8 approaching corresponding 
values for other baryons. Thus, not only the value but also the sign of the scalar self-energy 
for the S hyperon is sensitive to the strange quark content y. 




FIG. 3. Baryon effective masses as functions of y at the saturation density. 
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FIG. 4. The effective masses (a) and vector self-energies (b) for p, n, S+, S , H*^ and H at /3 = 1 
as functions of the baryon density. The calculation was performed at k = 8 and y = 0.08. 

F. Asymmetric nuclear matter 

Before discussing the asymmetric nuclear matter effects, let us recall general isospin 
symmetry relations which are automatically satisfied in the sum rule approach: 



sr(/3) 



/3), 
S^o(/3) =S|_(-/3), 



S^+(/3) 






m;(/3)=m:(-/3), 
m£4/3)=m*s_(-/3), 
m|„(/3)=m^-(-/3), 



(82) 
(83) 
(84) 



while A and S° self-energies do not depend on the isospin asymmetry. In the following, we 
will consider the extreme case of neutron matter with the asymmetry parameter (3 = 1 for 
p, n, S+, S^, H° and H~, while the effects of another extreme case of the proton matter 



with (3 = —1 can be obtained from the isospin symmetry relations (82-84). The effective 



masses and vector self-energies for p, n, S+, S , S° and S at /3 = 1 are shown in Fig. [4] as 
functions of the baryon density. 



The splitting of the effective masses for nucleons is determined by the I structure (78) 
resulting in negative proton-neutron mass difference Am* = m* — m* < in the neutron 
matter. This result is in agreement with leading OPE calculations in [IH]. Note, however. 



that inclusion of higher-order OPE contributions would provide Am*^ > which is expected 
in the relativistic approaches |19]. In contrast, the dominant contribution to the effective 



mass splitting for S and H hyperons comes from the q structure (77) resulting in relations 



m^_ > m^o > ""7.^+ and m|_ > m^o in the neutron matter. Numerically, the splitting at 
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the saturation density is of about 6%, 7% and 4% of the corresponding mass for A^, E and 
H baryons. 

The sphtting of the vector self-energies is dominated by the contribution of the p struc- 



ture (79). The relative strength of the splitting for different baryons mainly comes from the 
product p^^IsB, thus A^, E and S splittings should approximately scale as 6 : 7 : 1 in con- 
trast to the naive quark counting model which predicts the relation 1:2:1. Numerically, 
the vector self-energy splitting at the saturation density is of about 170, 260 and 40 MeV 
for A^, S and S baryons respectively in agreement with the approximate relation 6:7:1. 

IV. DISCUSSION 

In this section, we compare the obtained results with experimental data and alternative 
theoretical models. The nucleon in-medium properties have been already studied in detail 
in psl H9] . therefore we will concentrate on the hyperon case. We will discuss the solutions 
of the sum rule equations obtained in the approximation of the vacuum continuum thresh- 



olds (76) keeping in mind that exact solutions based on x fit (13) may lead to somewhat 
different results. 

Let us start with the A hyperon. The obtained scalar self-energy S^(po) ~ —210 MeV 
is in good agreement with the BHF calculations [T3] where the value m\/m\ = 0.84 was 
reported. Besides, the ChPT approach [20] provided the scalar self-energy about 55 MeV 
at 0.4po, which is close to the value 65 MeV obtained with the sum rule approach. As 
for the vector self-energy, the value S]((po) = 180 MeV was obtained providing a non- 
relativistic potential U\{po) ~ —30 MeV which is in surprisingly perfect agreement with the 
hypernuclear data [I] and the BHF calculations [15]. Recall, however, that the A effective 
mass is highly sensitive to k = {N\uu + dd\N) and hence to the value of the a.,^]si term. 
Changing the value from aT^^ = 45 MeV to 60 MeV would give much higher scalar self- 
energy and the potential would appear to be much deeper. 

On the other hand, the value of the effective mass m^ dramatically depends on the 
strange quark content y. Assuming y < 0.12, the value of m'^ would grow with density in 
agreement with the ChPT predictions [20]. However, for the larger strange quark content 
values, one would get the S effective mass decreasing with density in support of the BHF 
approach [15]. In any case, account of the large vector self-energy S^(po) ~ 0.3ms results in 
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the repulsive E hyperon potential at the saturation density in agreement with experimental 
data [2H3]. Touching the S isospin triplet in the neutron matter, we obtain the ordering 
S"*", S°, S~ in increasing mass shift similar to BHF calculations [SS]. Note, also, that the 
ChPT calculations [20j for the neutron matter provided anomalous ordering of the effective 
masses m'^o > ""^s- ^^^ ^° large isospin-symmetry violation effects. 

As for the S hyperon, the density dependence of the scalar self-energy appears to be 
similar to the nucleon case in contrast to the expectations from the naive quark model 
E|/E^ = 1/3. On the other hand, the S vector self-energy accounts about 1/8 of the 
nucleon self-energy. The large scalar self-energy in combination with the small E^ provides 
the attractive potential f/(po) about —200 MeV which is an order of magnitude larger 
than the value —18 MeV extracted from the experimental data [6'-'8]. Unfortunately, the S 
hyperon is not considered in ChPT or BHF approaches due to lack of SN scattering data. 



Note, however, that the Gell-Mann-Okubo-like formulas (65-66) could serve as an SU{3) 



motivated way for the tuning of the S-hyperon in-medium potential. 

The hyperon in-medium properties were also studied in the RMF framework [131 [H] 
where meson-hyperon coupling constants were fixed to reproduce the hyperon in-medium 
potentials at saturation density. This approach suffers from large ambiguities since the 
isoscalar a and u meson couplings appear to be highly correlated while the isovector meson 
couplings remain unconstrained by the hypernuclear data. As it was discussed in |48j, 
the lowest-order OPE terms, considered in this paper, correspond to exchanges by localized 
quark-antiquark pairs or effective vector and scalar mesons between baryons and nucleons in 
nuclear matter. Therefore one could use the obtained results as an input for the calculation 
of the effective meson-baryon coupling constants in the RMF approach. 

However, higher-order OPE terms appear to be numerically important. For example, 
inclusion of nonlocal vector condensate and four-quark contributions would subtract 60 MeV 
and 110 MeV from the lowest dimension value E]^(po) = 270 MeV ^48j. The scalar self-energy 
^Af(Po) = —140 MeV would remain almost unchanged because the four-quark condensates 
and nonlocal contributions would add about —100 MeV and 100 MeV, respectively |18] . 
Similarly, the higher-order OPE terms should play an important role in the hyperon sum 
rules. 
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V. SUMMARY 

The QCD sum rules provide a unique consistent formalism for investigation of the baryon 
octet in-medium properties. In contrast to other approaches, the QCD sum rules do not 
rely on phenomenological parameters of the baryon-meson interactions. 

In this paper, the baryon effective masses and vector self-energies were expressed in 
terms of a few in-medium QCD condensates of the lowest dimension which have been either 
calculated or related to the observables. It was shown, that the effective masses and vector 
self-energies in the baryon octet obey the relations similar to the Gell-Mann-Okubo mass 
formulas up to the linear S't/(3)-breaking terms. Moreover, the coefficients in the OPE 
terms provide a peculiar S'[/(3)-breaking pattern, e. g. vector self-energies in the symmetric 
nuclear matter are predicted to scale approximately as 1 : ^ : | : | for the A^, A, S and S 
baryons, respectively. 

Numerical studies for the A^, A, S and S baryon properties were carried out both in the 
symmetric and asymmetric nuclear matter in the approximation that in-medium effective 
continuum thresholds do not depend on density and remain equal to vacuum ones. The 
hyperon effective masses reveal a strong dependence on the values of a-j^N term and the 
strange quark content y which are known with poor accuracy. Nevertheless, the obtained 
effective masses and vector self-energies are in reasonable agreement with the results from 
other nuclear physics methods. 

The provided formalism can be extended to the case of matter composed of an arbitrary 
mixture of baryons, which is important in the calculations of the neutron star equation of 
state. Besides, contributions of the higher-dimensional condensates and radiative corrections 
could be included in order to improve the accuracy of the method. 
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